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ABSTRACT 

Ultraviolet (UV) galaxies have been selected from the GALEX deep imaging survey. 
The presence of a FUV-dropout in their spectral energy distributions proved to be 
a very complete (83.3 %) but not very efficient (21.4 %) tool for identifying Lyman 
Break Galaxies (LBGs) at z ~ 1. In this paper, we explore the physical properties of 
these galaxies and how they contribute to the total star formation rate. We divide the 
LEG sample into two sub-classes: red LBGs (RLBGs) detected at A = 24/im which 
are mainly Luminous IR Galaxies (LIRGs) and blue LBGs (BLBGs) undetected at 
A = down to the MIPS/GTO limiting flux density of 83 ^Jy. Two of the 

RLBGs are also detected at 70 fxm. The median SED of the RLBGs is similar (above 
A ~ 1/im) to that of a luminous dusty starburst at z ^ 1.44, HRIO. However, unlike 
local Luminous and Ultra Luminous IR Galaxies, RLBGs are UV bright objects. We 
suggest that these objects contain a large amount of dust but that some bare stellar 
populations are also directly visible. The median SED of the BLBGs is consistent with 
their containing the same stellar population as the RLBGs (i.e. a 250 - 500 Myrs old, 
exponentially decaying star formation history) but with a lower dust content. The 
luminosity function of our LBG sample at z ^ 1 is similar to the luminosity function 
of NUV-selected galaxies at the same redshift. The integrated luminosity densities of 
z ^ 1 LBGs and NUV-selected galaxies are very consistent. Making use of the RLBG 
sample, we show that star formation rates (SFRs) estimated from UV measurements 
and corrected using the IRX — (3 method provide average total SFRtot in agreement 
with the sum of the UV and infrared contributions: SFRjjv + SFRdust- However, 
IRX — /3-based SFRtot shows a large dispersion. Summing up the detected UV 
(1150A rest-frame) and IR-based star formation rates of the detected objects, we find 
that only one third of the total (i.e. UV -I- dust) LBG SFR resides in BLBGs and two 
thirds in RLBGs, even though most LBGs at z ~ 1 are BLBGs. On the other hand, 
the total SFR of LBGs accounts for only 11% of the total SFR at z 1. Finally, we 
observe a regular decrease of the luminosity ratio LdustI Lpuv from z = to z « 2 
for UV-selected samples. 

Key words: galaxies : starburst - ultraviolet : galaxies - infrared : galaxies - galaxies 
: extinction 



1 INTRODUCTION 

The evolution of the star formation density with redshift 
contains fundamental information from the perspectives 
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both of galaxy evolution and of cosmology. At intermediate 
redshift, star formation rates (SFRs) are estimated through 
a combination of UV, mid-IR, radio, and emission line tech- 
niques. Lyman Break Galaxies (LBGs) are used to estimate 
SFRs at high redshift. They are assumed to be related to 
FUV-selected samples selected via a color-color plane (e.g. 
Bouwens et al. 2006, Reddy et al. 2006). A major issue in 
determining the star formation density is to cross calibrate 
the several tracers used to estimate the SFRs. This paper 
tests whether the star formation rate estimated from the 
far ultraviolet (FUV) continuum (UV SFR or SFRj/y) of 
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UV-sclcctcd galaxies is equivalent to that from the infrared 
luminosity from dust estimated by integrating over the wave- 
length range 1 — lOOO/um (dust SFR or SFRdust, equivalently 
SFRT/i?). 

Even though we know that some dust is present in some 
LBGs (e.g. Chapman et al. 2002), it is difficult to get a di- 
rect estimate of how much in a statistically representative 
way because only a few LBGs at 2 > 2 have been detected 
in the FIR and in the sub millimeter (submm) ranges (e.g. 
Chapman et al. 2000). Even Spitzer has not detected a large 
sample of z > 2 LBGs in the regime where the dust emis- 
sion is the main component of the spectrum (e.g. Huang et 
al. 2005). Adelberger & Steidel (2000) concluded that the 
bulk of the submm background is produced by moderately 
obscured galaxies similar to the ones already detected in UV- 
selected surveys. Most of the dust-to-FUV luminosity ratios 
{Ldust/Lpuv) for these galaxies are estimated from the ob- 
served optical and by using the UV slope vs. Ldust/Lpuv 
(IRX-/3) relation defined by Meurer, lieckman & Calzetti 
(1999). However, the validity of this law is debated (Bell 
2002, Buat et al. 2005, Goldader et al. 2002). Nevertheless, 
using this law, dust luminosities are often estimated with- 
out FIR data (e.g. Adelberger & Steidel. 2000). Can we trust 
these results? Sawicki (2001) concluded that "it is unlikely 
that a consensus about the amount of stellar light inter- 
cepted by dust will be reached on the basis of rest-frame 
UV and optical data alone" . 

We therefore need rest-frame FIR data to advance our 
knowledge of the dust content of UV-selected galaxies. In 
this paper we approach this goal by observing LBGs at 
2: ~ 1 where we can directly evaluate the dust luminosity 
from 24^m data. Burgarella et al. (2006b) studied an ini- 
tial sample of LBGs in the redshift range 0.9 2 ^ 1.3 
identified from GALEX data in the ultraviolet; a portion 
of this sample was also detected with Spitzer /MIPS in the 
IR. The remainder are below the limiting MIPS flux density 
of 83fiJy. Although this is a true Lyman Break sample at 
2 ~ 1, the selection is dift'crcnt from LBGs at z ^ 3 be- 
cause the effect of the intergalactic medium is much lower 
at « 1 than at higher redshift. The same phenomenon 
applies, although at a lower level, for BM/BX galaxies se- 
lected at 2 ~ 2 (e.g. Reddy et al. 2006). Thus, although exact 
analogs of the high redshift LBGs cannot be used, the main 
objective in the previous work was to define a rest-frame 
FUV-selected sample with a measured estimate of the dust 
luminosity from IR data. 

In this paper, we expand upon this approach. We de- 
scribe the selection of the sample in Section 2. Section 3 is 
devoted to the building and the analysis of the spectral en- 
ergy distributions (SEDs) of LBGs (i.e. red LBGs (RLBGs) 
and blue LBGs (BLBGs)) and to their comparison to the 
SEDs of LBGs at « « 3. Then, in Section 4, we will derive 
the luminosity function of the z ^ 1 LBG sample and com- 
pare it to the luminosity function of an ultraviolet-selected 
sample. In Section 5, we will address the question of the to- 
tal star formation density at 2; ^ 1 and whether this total 
star formation density can be evaluated from an ultraviolet- 
selected sample, an infrared-selected sample or whether both 
of them are required. In Section 6, we will study the evolu- 
tion of the LcLust/ Lfuv ratio for UV-selected samples or, in 
other words, of the evolution of the dust attenuation with 
redshift. We will finish this paper in Section 7 by estimat- 



ing how we can detect LBGs in the infrared with facilities 
planned for the next decade. 

We assume a cosmology with Hq = 70 km.s~^ .Mpc~^ , 
Qm = 0.3 and ilvAC = 0.7 in this paper. 



2 DEFINITION OF THE UV AND LBG 
SAMPLES 

Our first paper (Burgarella et al. 2006b) on the LBG sam- 
ple at 2 ~ 1 used GALEX pipeline data. Here, we use the 
deeper (76444 sec) NUV image available in GALEX Release 
2 (GR2) whose central coordinates are a = 03h32m30.7s, 
and 5 = -27deg52'16.9" (J2000.0). An analysis of GALEX 
deep images shows that galaxies overlap, producing an ap- 
parent brightening of some of the detected objects in the 
tables available from the Multi Mission Archive at STScI 
(MAST). This effect starts to apply to objects brighter than 
NUV « 22.0. 

To deal with these effects systematically, we put all of 
our images at a common resolution and extracted photom- 
etry by means of point spread function fitting. The spatial 
resolution of the GALEX images was slightly degraded by 
applying a median 3x3 filtering to the FUV and NUV im- 
ages. The GALEX "beam" is taken to be the solid angle 
of the la radius circle of a Gaussian point-spread function 
(PSF) : Qbeam = ttct'^ as defined in Hogg (2001). Note that 
for our GALEX analysis, Ofwhm = 6 arcsec which trans- 
lates into cr = OpwHM l^.'ih — 2.55 arcsec. The final angular 
resolution in both bands therefore, matches the MIPS 24/im 
angular resolution. 

Confusion noise is significant in these images (see Dole 
et al. 2004 for 24^m). The size of the filtered GALEX beam 
Hbeam = 4.81 X 10""' Sr. Slucc the total field of view of 
GALEX corresponds to VIgalex ~ 3.45 x lO""* sr, we obtain 
a confusion limit at 3(j: 

(s/6)co„y = — — — '^"""'''^'^ = 0.063 source per beam 

or ~ 16 beams per source at the 80% completeness limit 
defined below for the NUV-selected sample. This is deeper 
than the criterion published by Hogg (2001) of 20 beams 
per source but shallower than the ~ 12 beams per source 
estimated by Dole et al. (2004) and Jeong et al. (2006). 
Takeuchi & Ishii (2004) have published a relation to estimate 
the confusion limit at 5cr. Our number counts can be fitted 
by a power law with a slope 7 ~ 1.45. Using e = 2 and 
a = 5 as recommended by Takeuchi & Ishii (2004), we can 
compute: 

{s/b)conf = ^^22"^^ ~ 0.044 source per beam 

or 23 beams per source. Accounting for the difference in 
detection level {3a for the former and 5a" for the latter), the 
two values are consistent. 

We used DAOPHOT for photometry in these images 
with a PSF built from 10 bright point sources in the FUV 
and NUV. DAOPHOT was developed for point sources. Le 
Floc'h et al. (2005) show it is suitable to measure 24/jm 
flux densities of Spitzer/MIPS sources. Our UV sources are 
restricted to the redshift range 0.9 < 2; ^ 1.3. De Mello et 
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al. (2006) showed that UV-selected galaxies in the redshift 
range 0.8 s: z ^ 1.3 have = 1.6 ±0.4 kpc i.e. FWHM « 6 
kpc. At 2 ~ 1 this size corresponds to 0.7 arcsec. This is weU 
below the angular resolution of our images, so we can safely 
use DAOPHOT for flux estimation in the targeted redshift 
range. 

We detect about 10^ objects in the NUV down to 3 
a with the DAOPHOT DAOFIND task, which is in agree- 
ment with number counts by Xu et al. (2005). Using the 
DAOPHOT ADDSTAR task, we estimate the 80% com- 
pleteness limit to be NUV=26.2 for the NUV-selected sam- 
ple. The NUV-selected sample used hereafter corresponds 
to all objects with a NUV magnitude NUV ^ 26.2; it con- 
tains 45366 objects in the GALEX fleld of view but only 
part of them have a COMBO 17 counterpart (i.e. a red- 
shift) mainly because most GALEX objects lie outside the 
COMBO 17 field of view. This source density corresponds 
to ~ 16 beams per source, that is, just at the 3-a confusion 
limit. 

To check the validity of our approach, we also apply 
DAOPHOT PHOT aperture photometry (integrated over 3 
arcsec with an aperture correction estimated from the PSF) 
to a set of isolated objects (i.e. no COMBO 17 counterparts 
within a 4-arcsec radius) that we compare to DAOPHOT 
magnitudes for the same objects. The difference in magni- 
tude (DAOPHOT - PHOT) amounts to 0.17 ± 0.26 which 
is about the uncertainty of a typical GALEX measurement 
(Morrissey et al. 2005) given that we reach very faint magni- 
tudes, confirming that DAOPHOT provides good estimates 
for the GALEX magnitudes. 

Next, we measure the FUV flux of those NUV-selected 
objects by applying DAOPHOT on the FUV images at 
the location of NUV-selected objects. Tests show that our 
FUV data are 80 % complete down to FUV=26.8. Then, 
we perform a cross-correlation with COMBO 17 (Wolf et al. 
2004) within a radius r=2 arcsecs. We find a counterpart in 
COMBO 17 for about 71% of the NUV-selected sample. Ta- 
ble 1 provides the calibrated output from DAOPHOT (FUV 
and NUV) and the aperture photometry within a 3-arcsec 
radius for the resulting 7309 objects inside the 0.263 deg^ 
m common with COMBO 17. H Some of the sources have 
several possible counterparts (several NUV counterparts for 
one COMBO 17 source or vice- versa). An object with an 
observed color FUV — NUV ^ 2.0 is classified as a member 
of our LBG sample if the redshift of the counterparts (some- 
times several within 2 arcsecs) from COMBO 17 is in the 
range 0.9 ^ 2 ^ 1.3 (to be safe, we discard objects with vari- 
ability flags larger than 8 and objects classified as "Star", 
"WD", "QSO"; "QSO (gal?)" and "Strange Objects"). At 
this point, we have 420 LBGs (2 with four COMBO 17 
counterparts, 40 with two counterparts and 378 with one 
counterpart). We estimated that the total completeness of 
the NUV-selected sample is about 80%. Correcting for this 
80% UV completeness, the observed surface density of NUV- 
selected galaxies (including LBGs) down to NUV = 26.2 in 
the 0.26 sq. deg. is 6850 deg~'^. 



^ Some additional information such as the redshifts and the R- 
magnitudes extracted from COMBO 17 are also provided. How- 
ever, the original COMBO 17 database should be consulted for 
complete information, using the listed COMBO 17 identification. 



Table 2. The efficiency of our selection to identify 0.9 ^ ^ ^ 1.3 
galaxies through the detection of FUV dropout galaxies is 21.8 
% in the range 20.8 NUV ^ 26.2 but the completeness is very 
high at 83.3 %. 



^^UV-set,allz 


-^break ,allz 


NLBG,0.f)^z^l.3 


Efficiency 


8206 


2170 


474 


21.8 % 


Nuv-sel,0.9^z^l.3 






Completeness 


1943 


568 


473 


83.3 % 



The 80 % completeness is reached at FUV=26.8. Even 
though we still measure fluxes below this limiting magni- 
tude, LBGs are missed for galaxies with NUV ^ 24.8 be- 
cause FUV-NUV cannot be measured. However, if we re- 
strict ourselves to NUV < 24.8 where we are able to esti- 
mate FUV — NUV without the above restriction i.e. about 
80% completeness, the respective surface densities of NUV- 
selected objects and LBGs are 1766 deg~^ and 1320 deg^^ 
i.e. about 75% of the NUV-selection are LBGs. A deeper 
analysis is postponed to Sect. 4. 

Table 2 shows the completeness and the efficiency of 
the criterion FUV — NUV ^ 2 for identifying star-forming 
galaxies in the redshift range 0.9 ^5 z ^ 1.3. This method 
is very complete whenever the FUV — NUV color is fully 
measurable i.e. when NUV < 24.8 as described above. A 
more efficient method over all redshifts would need to use 
color-color diagrams as suggested by the statistics presented 
by Adelberger et al. (2004). However, this is not the goal 
of this paper, where COMBO 17 redshifts are available as 
inputs to the sample selection. 

Finally, we carried out cross-correlations to complete 
the wavelength coverage with the UBVRIJK ESO Imaging 
Survey (EIS), four IRAC bands and finally with MIPS at 24 
Hm and 70/xm. We stress that no additional selection effects 
enter into the final LBG sample since all the objects are kept 
in the following analysis with their own detections / non de- 
tections depending on the band. We have therefore, at the 
end of the cross-correlation process, the same number of UV 
objects as after the UV selection. Table 3 lists the best red- 
shift as in Table 1, the UV 1800 A luminosity estimated by 
interpolating the observed NUV and U-band u and the 
total IR luminosity evaluated from the 24/im flux density 
and Chary & Elbaz (2001) calibration. An interesting point 
is that 62 of the 420 LBGs (i.e. 15%, but 17% if we restrict 
our sample to NUV < 24.8 as above) have a 24^m MIPS 
counterpart down to the MIPS/GTO limiting flux density 
of 83 iJ,Jy (all of them except one within r=2 arcsecs); we 
will call them Red LBGs (RLBGs hereafter). Huang et al. 
(2005) observed LBGs at z ~ 3 with Spitzer and, excluding 
objects with possible contamination from AGN, they also 
classified about 15% as IR-luminous LBGs. Only two RL- 
BGs are detected at 70fim but none of them are detected 
at 160/im. The remaining LBGs, undetected in the MIPS 
image are called Blue LBGs (BLBGs). Note that this is a 
practical notation and it does not imply a physical differ- 
ence so far. The observed surface densities are 1657 deg~^ 
for the much more numerous BLBGs, and 339 deg~'^ for the 
RLBGs down to NUV = 26.2. 
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Table 1. Ultraviolet-selected sample down to NUV=26.2 within COMBO 17 (C17) field of view with magnitudes estimated from 
DAOPHOT. We stress that the magnitudes of resolved objects must be used with care because DAOPHOT is designed for point sources. 
However, for galaxies beyond z ^ 0.5, DAOPHOT photometry is valid (see text for a more detailed discussion). 46321 rows for 45366 
iduv are listed out of which 38057 correspond to GALEX objects without COMBO 17 counterparts, 6404 GALEX objects have a single 
C17 counterpart and a few ones with up to 4 counterparts (859 X 2, 42 X 2 and 4 X 4). Several rows with the same idyy are listed 
whenever several C17 counterparts are found. For each row, columns give (1) the number of C17 counterparts (1), the UV id (2) and 
coordinate (3,4), the DAOPHOT FUV and NUV ABmag from PSF fitting with their estimated uncertainties (5-8); the C17 id (9), the R 
magnitude and uncertainty (10,11), the photometric redshift from C17 with a quality flag (0 means from MC.z-ml and (1) from MC-Z 
in C17) and finally the magnitude estimated from PHOT (3 arcsec aperture). 



ncl7 


iduv 


RAc/v 




FUV 


errpuv 


NUV 




idci7 






z 




NUVpft„t3 


errjvc/Vp^, 


1 


41129 


52.946268 


-27.968764 


23.967 


0.050 


23.636 


0.038 


12685 


21.311 


0.011 


0.187 


1 


23.570 


0.027 


1 


41130 


52.934597 


-27.968940 


27.583 


0.221 


25.409 


0.039 


12616 


21.827 


0.024 


0.729 


1 


25.445 


0.144 


2 


41133 


52.852864 


-27.968810 


26.600 


0.111 


24.159 


0.038 


12556 


22.492 


0.021 


0.616 


1 


24.164 


0.049 


2 


41133 


52.852864 


-27.968810 


26.600 


0.111 


24.159 


0.038 


12627 


21.465 


0.017 


0.685 


1 


24.164 


0.049 


1 


41134 


52.830309 


-27.968789 


24.451 


0.060 


22.685 


0.046 


12683 


19.897 


0.004 


0.535 


1 


22.707 


0.016 



Table 3. The LBG sample at 0.9 ^ 2 ^ 1.3, as defined in the text, is listed here (466 rows). Rcdshifts arc extracted from COMBO 
17 and the redshift quality flag Zq^al corresponds to '1' if z is extracted from column "MC_z" and '0' if z is extracted from column 
"MC_z_ml" in Wolf et al. (2004). No values '-99' are listed in the columns if the LBG is not detected in the U band (i.e. Lisoo) or at 
24/im (l-iduat) down to 83 fJ.Jy, i.e. BLBGs. By deflnition, all objects are observed in NUV and should have a Ln.-jo value. Objects with 
^ ^dust value are RLBGs. Aigoo is computed from the l-idust I Lieoo ratio (Lieoo interpolated from the observed NUV and U bands). 



idf r 


«Arv 


D.vrr 


L'iir.i) 




L,/„., 


z 


z,,„„; 


Aieoo 


46333 


53.106075 


-27.918961 


9.800 


10.193 


11.249 


0.900 


1 


2.092 


39913 


53.252274 


-27.980711 


9.540 


10.006 


-99.000 


0.900 





-99.000 


77655 


52.931938 


-27.613210 


10.054 


10.185 


11.008 


0.903 


1 


1.587 


50879 


52.957910 


-27.876471 


9.580 


9.482 


-99.000 


0.904 


1 


-99.000 


53306 


53.220624 


-27.853428 


10.190 


10.084 


-99.000 


0.905 


1 


-99.000 



3 THE SPECTRAL ENERGY DISTRIBUTIONS 
OF LYMAN BREAK GALAXIES 

In this section, we analyse the median SEDs of RLBGs and 
BLBGs. In a previous paper (Burgarella et al. 2006b), we 
had stacked the 24 /im images of BLBGs in a first z ~ 1 
LBG sample. Although we have improved our UV photo- 
metric measurements by using DAOPHOT, we consider that 
the stacked 24/xm flux densities for RLBGs given in Bur- 
garella et al. (2006b) are valid because i) most of the LBGs 
in the first sample are included in the present one and ii) 
the stacking process is designed to reduce the influence of 
objects around the central detection (i.e. the LBGs) to en- 
hance the signal-to-noise ratio of individually undetected 
BLBGs. The effect of the confusion is therefore reduced. 
Burgarella et al. (2006b) had estimated an average flux den- 
sity of IS/uJj/ at the observed wavelength of 24 /jm for those 
BLBGs. Table 4 gives the absolute B magnitudes Mb (from 
the observed I band) of BLBGs and RLBGs and their FUV 
absolute magnitudes. As order-of-magnitude comparisons, 
our LBGs are more luminous than Blue Compact Galax- 
ies (e.g. Mb < -18.5 fi:om Noeske et al. 2006); BLBGs are 
fainter than and RLBGs similar to LBGs at « ~ 3 (e.g. 
< MiTooiBLBG) >= -21.0 ± 1.0 in Sawicki & Thompson 
2006). 

Ideally, to relate our galaxies to each other it would be 
better to use the redshift of each galaxy to calculate K cor- 
rections and put them all in the same frame before building 
the template SED. However, since we cannot simply inter- 
polate between observed measurements because of high con- 



Table 4. Average absolute magnitudes of the 2 ~ 1 LBG sample 





Mb 






Blue LBGs 


-21.0 ± 


1.0 


-20.3 ± 1.0 


Red LBGs 


-22.4 ± 


1.0 


-20.9 ± 1.1 


AU LBGs 


-21.0 ± 


1.1 


-20.3 ± 1.1 



trast spectral features (for instance an interpolation between 
the observed 8 nm and the observed 24 /irn bands would 
produce catastrophic results), we would need to use some 
fiducial SED for interpolation in redshift. Since there axe 
no observed FUV-to-FIR SEDs built from a large sample of 
high redshift galaxies, we would have to rely on models. We 
prefer to avoid the resulting unknown uncertainties. Given 
our small redshift range, K-corrections should be relatively 
small. We plot in Figs. 1-3 the extent of the wavelength 
range (except in Fig. 3 hidden by symbols) due to the range 
in redshifts, showing that the errors are small from assuming 
an average redshift of z = 1.1 for all galaxies. 

3.1 Red Lyman Break Galaxies (RLBGs) 

We can study for the first time a large sample of LBGs in 
the wavelength range where the emission by dust is pre- 
dominant. Moreover, two IR-bright LBGs in our sample 
have a significant piece of SED in the FIR with data points 
at 24/im and 70fj,m. Interestingly enough, these two ob- 
jects are consistent with being simultaneously ULIRGs (i.e. 
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LiR > W^^Lq) and very blue (UV slope f3 2) UV lu- 
minous galaxies (UVLG) with moderate dust attenuation 
(Afuv ~ 1.8). They will be described in a forthcoming 
paper (Burgarella et al. 2007). The other RLBGs are not 
detected at 70/im. 

To build the UV-to-mm SED of RLBGs presented in 
Fig. 1, we proceed as follows: 

• We selected the RLBGs that are detected at all wave- 
lengths from the observed NUV to 24fj,m (i.e. 37 objects) . If 
we had computed the median SED of the full LBG sample, 
we would have derived a SED biased toward the brightest 
UV objects (and left undefined how to treat the many up- 
per limits at 24//m) . By forcing the sample to have measured 
data points over the full wavelength range, we moderately 
bias the resulting SED toward red objects, which is consis- 
tent with our goal to characterize red LBGs. 

• We normalize the median SED of two RLBGs detected 
up to 70/Lim to the median SED of RLBGs in the observed 
IRAC 8/im band. This normalization demonstrates that the 
24pm flux densities of these two RLBGs are very consistent 
with the 24/xm flux density of the median SED of the RL- 
BGs. The values at other wavelengths are fairly consistent 
as well. The two SEDs are very likely extracted from the 
same parent RLBG population. Therefore, we make the as- 
sumption that the normalized 70jj,m flux density of the two 
RLBGs detected at this wavelength can be used to extend 
the template SED of RLBGs. However, we must remember 
that the scatter in the two 70/xm data points is very likely 
not representative of a larger sample. 

• We normalized the SED of HR 10, a dusty and luminous 
{Ltir ~ 9 X lO^^L©) starburst at z=1.44 (Stern et al. 2006) 
to the observed IRAC 3.6/^m band where the output should 
be dominated by stellar photospheres {Xrest ~ 1.8/im). It 
is interesting that the SED matches very well that of the 
median RLBG SED for rest-frame wavelengths in the range 
1.8/ttm - 70fj,m (Fig. 1). Again, it seems that we are deal- 
ing with the same kind of objects (e.g. Elbaz et al. 2002). 
Therefore, we tentatively use the SED of HR 10 to extend 
our median SED to wavelengths longer than 70/ttm up to 
1mm. 

The resulting RLBG template is listed in Table 5 and plot- 
ted in Fig. 1. The RLBG SED does not exhibit a power-law 
shaped continuum characteristic of AGN-dominated emis- 
sion. Indeed, Ivison et al. (2004) showed that strong AGNs 
have nearly a constant slope from 2 to 10 iJ,m while a distinct 
minimum in the 3 to 4 fj,m range is observed for staxburst- 
dominated galaxies, and this property is being used to iden- 
tify AGN through IRAC colors (e.g.. Lacy et al. 2004; Stern 
et al. 2005; Alonso-Herrero et al. 2006). A key diagnostic 
for distinguishing AGN-dominated and starburst-dominated 
galaxies is a change of slope between the stellar and dust 
continua at 3-6 /xm, which is one of the most noticeable 
features in the SEDs of our RLBGs. 

However, even though HR 10 bears strong similarities 
in the IR to our RLBG population, its SED is very faint in 
the optical and UV and we cannot simply state that RL- 
BGs are HR 10 - like objects. The RLBGs might represent 
a type of galaxy at 2 ~ 1 that could be considered to be 
a missing link between the two high redshift populations 
of galaxies: the blue and almost dust-free LBGs and the 
IR/sub-mm bright LIRGs/ULIRGs. We might be observing 



complex objects where one component is emitting in the UV 
and the other in the IR. A related object might be VV114 
(Le Floc'h et al. 2002), which is also considered to be a local 
counterpart of distant LBGs. Do we see different classes of 
unrelated objects or is there an evolution from one class to 
the other where RLBGs represent a temporary phase? We 
do not have enough information, yet, to decide and a mor- 
phological analysis of the full sample is under way from HST 
imaging to provide clues toward or against this hypothesis. 

To understand this median SED, we model it using PE- 
GASE (Fioc & Rocca-Volmerange 1997) in the UV + Op- 
tical + NIR (UVONIR) and several modified blackbodies 
in the FIR and sub-mm ranges. In the UVONIR,, the best 
match is for a 500 Myr-old model (exponentially decaying 
SFR *(t) = exp{-t/r) with r = 7 Gyr, 1/10 solar metallic- 
ity and a Salpeter IMF). To redden the model, we assume 
an attenuation law with Afuv oc (X/Xfuv)" where Afuv 
is the FUV dust attcmiation. Since all RLBGs have Afuv 
estimated from IR data, we used their mean Afuv = 2.45 as 
a constraint. The slope, n, is a free parameter and n = —0.7 
provides the best match to the observed SED. The UVONIR 
SED model is normalised to the observed R-band photomet- 
ric point. 

In the MIR/FIR, we use four modified blackbodies with 
an emissivity index set to e = 1.5 (Klaas et al. 1997). These 
blackbodies are not meant to reproduce the detailed emis- 
sion spectrum in the PAH region (i.e. below about 12/im, 
e.g. Marcillac et al. 2006, Smith et al. 2006) but orfly the 
underlying continuum. As a consequence, dust luminosi- 
ties estimated from this approach are, strictly speaking, 
only lower limits. The temperatures of the four blackbod- 
ies are T = 300 K, 120 K, 45 K and 20 K. The normal- 
ization of the hottest one is fixed by setting the sum of 
the UVONIR model plus the 300 K model to the observed 
IRAC point at 8.0 nm. The warm blackbody is constrained 
by the 24/im MIPS point. The cool blackbody is fixed by 
the observation of the two RLBGs detected at 70^m and 
finally the cold blackbody by the 350 /irn measurement 
of HR 10 (see above). The total IR luminosity from 1 to 
1000 iim amounts to Log{LTiR) = 11.5 from the median 
RLBG SED. These RLBGs are therefore LIRGs on aver- 
age. The UV luminosity integrated from 0.1 pm to 1.0 fim 
amounts to Lfuv = 10.9 but to estimate the Ltir/ Lfuv 
ratio, we estimate A. /a (A = I8OOA) = 10.2. These objects 
are therefore not UV luminous galaxies (UVLGs) on aver- 
age. The resulting value of Loq^Ltir)/ Lfuv = 1.3 converts 
into Afuv = 2.5 using the calibration of Burgarella et al. 
(2006a). This value is in excellent agreement with the dust 
attenuation estimated from the fit. 

3.2 Blue Lyman Break Galaxies (BLBGs) 

To build the UV-to-mm SED of the BLBGs presented in 
Fig. 2, we proceed as follows: 

• The median SED of BLBGs is computed for BLBGs 
detected at all wavelengths from the observed NUV to the 
IRAC 8/xm band (i.e. 68 objects). 

• The 24/im point is the result from the stacking carried 
out in Burgarella et al. (2006b). 

• The median SED of RLBGs presented above is super- 
imposed onto the median SED of BLBGs after normalizing 
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it to the observed IRAC 8/xm band, where stellar emission is 
still predominent. The match is fairly good from the IRAC 
3.6/um band to 24fj,m and we use this normalized RLBG SED 
to extend the BLBG one, which means that we make the as- 
sumption that the temperatures of dust grains in BLBGs arc 
the same as in RLBGs. Only the amount of dust attenuation 
differs. 

As for RLBGs, the median SED of the BLBGs does not 
show a power-law shaped continuum characteristic of AGN- 
dominated emission in the 3-4 /jm range. As we showed 
for the RLBGs, the BLBGs should be starburst-dominated 
galaxies. 

The UVONIR part of the SED is fitted again by a model 

from Fioc & Rocca-Volmerangc (1997). The best match is 
for a 250-Myr model with the same global characteristics as 
the RLBGs except for the amount of dust attenuation. For 
the BLBGs, we do not have any individual measurements 
of Ldust- A constraint can be set by the upper limit at 83 
fijy flux density at 24 ^m, which can be converted into an 
upper limit of Afuv ~ 2.1. The best matching model cor- 
responds to Afuv ~ 1.8 with the same slope n = —0.7 as 
for the RLBGs. We conclude that the stellar populations 
in the BLBGs are only slightly younger than in the RL- 
BGs and that the main difference between the two types 
comes from the amount of dust attormation. The total IR 
luminosity of the median BLBG SED from 1 to 1000 fim 
amounts to Loq^Ltir) = 10.9. A median BLBG is there- 
fore not a LIRG. The UV luminosity integrated from 0.1 
/im to 1.0 /im amounts to Lfuv ~ 10.6. Again, to esti- 
mate Ldust/ Lfuv , we estimate A f\{X = I8OOA) = 10.1, 
which means that BLBGs are not UVLGs. The resulting 
value of LogLdust/Lpuv = 0.9 converts to Afuv = 1-8 
which, again, is consistent with the above value and with 
Schiminovich et al. (2005). 

3.3 Comparison with LBGs at a ~ 3 

Our two SEDs for BLBGs and RLBGs axe compared in Fig. 
3 to an average SED estimated for LBGs at a 3 in Forster- 
Schreiber et al. (2004) and normalised at 2750 A. The grey 
cloud represents the standard deviation in Forster-Schreiber 
et al.'s LBG sample. Although marginally consistent with 
the latter, the SED of the RLBGs is too red in the UV 
and has an excess in the red part of the SED. Older stellar 
populations and/or more likely larger dust attenuation pre- 
sumably are the origin of this difference. On the other hand, 
the SED of BLBGs has strong similarities to « ~ 3 LBGs: 
all the data points appear to be consistent, within the un- 
certainty, with the hypothesis that the two populations are 
extracted from the same sample. More data are needed in 
the rest-frame NIR and in the MIR/FIR to provide infor- 
mation on the evolved stellar population and on the amount 
of dust to reach a definite conclusion. However, the present 
data and the results of Huang et al. (2005) suggest that we 
might observe LBGs at a ~ 3 with the same SED char- 
acteristics (shape and distribution) as LBGs in the 2; ~ 1 
universe. Moreover, the LBG selection seems to favour star- 
bursting galaxies with low dust attenuation whatever the 
redshift, that is BLBGs as defined in our sample. Neverthe- 
less dustier LBGs (i.e. RLBGs) exist at ^; ~ 1 and z ~ 3 
with a higher total SFR than BLBGs. However, we saw in 



Table 4 that LBGs at z ~ 3 are brighter than BLBGs but 
have magnitudes similar to RLBGs. This means that a sim- 
ple identification of BLBGs with high redshift LBGs is not 
strictly valid. 



4 THE ULTRAVIOLET LUMINOSITY 

FUNCTION OF LYMAN BREAK GALAXIES 

The LBG selection process preferentially detects stax- 
forming galaxies with low dust attenuation (see also, e.g., 
Steidel et al. 1996). This is generally true for UV-selected 
galaxies whereas IR-sclected galaxies tend to be star-forming 
galaxies with a high Ldust/ Lfuv ratio, that is, high dust 
attenuation (Iglesias-Paramo et al. 2006, Buat et al. 2006, 
Burgarella et al. 2006a). Xu et al. (2006) suggest that these 
differences are mainly caused by selection effects: both UV- 
sclected and IR-sclected galaxies would be extracted from 
the same population of star-forming galaxies. Buat et al. 
(2006) reach a somewhat similar conclusion. What can we 
say about high redshift LBGs? 

LBGs are selected through colours. Therefore, strictly 
speaking they do not form a genuine UV-selectcd sample. 
Their spectral break is caused by the absorption of UV pho- 
tons with A < 912A in the outer atmospheres of massive 
stars, in interstellar H I gas, and in intervening H I gas along 
the line of sight. However, in the case of star-forming galax- 
ies, the flux redward of the Lyman break (in the rest frame 
UV) is high, and high redshift LBGs axe often assumed to 
be similar to a UV selected sample (Giavalisco 2002). High 
redshift LBGs are therefore used to estimate the UV star 
formation density. However, the intergalactic medium has 
little effect in producing the Lyman break in the SED of 
LBGs at 2; ~ 1 and our objects are true Lyman break galax- 
ies for which the break is produced by material inside the 
observed galaxies. We need therefore to check whether our 
LBG selection is equivalent to a UV selection before we can 
put local and high redshift LBG studies on the same ground 
to compare them. 

To address these questions, we need to compare the lu- 
minosity function (LF) of a UV-selected sample to the LF of 
a LBG sample. We address this point at « ~ 1. Arnouts et al. 
(2005) computed the LF of UV-selected galaxies at 150nm 
and z ^ 1. We evaluate the LF of our NUV-selected sam- 
ple using the new DAOPHOT photometry and the 1/Vmax 
and C— methods, to which we applied a K-correction to get 
LFs in the same wavelength range as for the higher redshift 
studies (e.g. Madau et al. 1996, Steidel et al. 1999, Sawicki 
& Thompson 2006). Both LFs at 2: ~ 1 arc very consistent 
(Fig. 4). We also compute the Far-UV (FUV) LF of our LBG 
sample (using 1/Vmax and C— ). To select the z ~ 1 LBG 
sample, we use a color selection FUV — NUV > 2. This 
means that there is a secondary selection in addition to the 
primary FUV selection on the sample. This secondary selec- 
tion is properly addressed by the Kaplan-Meier estimator. 
By this method, we can utilize the information contained in 
the upper limits (see a more detailed description in Xu et al. 
2006), and correctly estimate the LBG luminosity function. 
We compare it to the two FUV LFs in Fig. 4. 

The UV luminosity density ppuvi^SOOA) from Schimi- 
novich etal. (2005) is p|^v(1800A) = 9.3xlO'';i[L©Mpc-^]. 
If we restrict Arnouts' LF to our sample range, we obtain 
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Table 5. Measured median SED of 37 RLBGs and 68 BLBGs u (erg.cm-^.s-l). For RLBGs, the 70 IMm data are extracted for two 
RLBGs only. The 24 fim data for BLBGs have been estimated by stacking BLBGs (Burgarella et al. 2006b). 
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Figure 1. The average SED of RLBGs is modelled using a reddened PEGASE model in the UVONIR, and modified blackbodies in the 
FIR and submm ranges. Good agreement with the SED of HRIO exists in the IR, showing that we see very dusty galaxies. However, 
unlike HRIO, RLBGs are also detected in the UV. 
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Figure 2. The average SEDs of BLBGs looks similar to that of RLBGs but we have to assume a lower dust content to explain the 
UVONIR SED. This is in agreement with the non direct detection at 24/im. The value at this wavelength is provided by stacking. 



\CZ^ LBGs at z~3, Forster-Schreiber et al. (2005) 
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Figure 3. The UV/optical SEDs of the two classes of LBGs are noted as in Figs. 1 and 2. They are compared to the average SED of 
LBGs at 2; ~ 3 by Forster-Schreiber et al. (2004) located in the grey area. The comparison suggests that the latter are more similar to 
BLBGs than to RLBGs on average. 
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Figure 4. The luminosity function of all LBGs estimated from 
the 1/Vmax method and confirmed by the C— method. At 
Log{Lpuv) > 2 X 10^" i.e. UVLGs, this luminosity function is 
almost the same as the NUV luminosity function of UV-selected 
galaxies at 2 = 1 from our own measurements and from Arnouts 
et al. (2005). 



sel LF gives 



5.1 X IQ'hlLQMpc-^] while our NUV- 



Pf(7v(1800A) = 4.0 ± 0.6 X lO^ft[L0Mpc~% 
The difference amounts to ~ 20%. For LBGs, 

PFa?(1800l) = 3.4 ± 0.7 X lO^/i[L0Mpc"^] 

which represents 85% of the total pF(7v(1800A) of NUV- 
selected galaxies in the considered luminosity range. The 
given uncertainties are 20%, so the two values are consistent. 
We can conclude that at z ~ 1, our LBG selection (which 
might be different from a higher redshift LBG selection i.e. 
without screening objects with high Apuv) is similar to a 
UV selection. 

Since most of our LBGs (83%) belong to the BLBG 
class, the difference between the LFs of BLBGs and LBGs is 
small. As expected, the LF of RLBGs (~17 % of the sample) 
falls well below. While the shape above Log{LFuv / Lq) = 
10.3 (i.e. UVLGs) of the RLBG LF is the similar to that of 
all LBGs, the RLBG LF seems to flatten at low luminosity. 
The contributions from blue and red LBGs to the LBG UV 
luminosity density are as follows: 

pfbv°{i800A) = 2.9 ± 0.9 X w''h[LQMpc~^] 
i.e. 72.5% of the FUV flux density, and 



Pfuv° i'^SOOA) = 2.7 ± 1.1 X 10^h[LQMpc-^] 
i.e. 7% of the FUV flux density. 



5 CONTRIBUTIONS OF UV- AND 
IR-DEDUCED SFR TO THE TOTAL STAR 
FORMATION DENSITY 

Assuming a 250-500 Myr exponentially decaying star forma- 
tion rate from PEGASE, in agreement with the SED fitting 
in the previous section, we can estimate SFRnso for the 
UV-selected sample and for LBGs: 

SFRiiso = i^L^ [erg s''^]/ 1.8 x 10"^° Mq yr"^ Mpc^ 

In determining this value, we used our 80 % complete NUV 
measurements, which apply to a rest-frame wavelength of 
1150A. This wavelength corresponds to the observed V-band 
at z 3 — 4. Since we will make comparisons within our own 
dataset, adopting this wavelength reduces the uncertainties 
due to K-corrections. 

Using the observed 24/xm flux density we can estimate 
the luminosity from emission by dust, Ltir, from the Chary 

6 Elbaz (2001) calibration. However, it is necessary to eval- 
uate how the aromatic bands could affect the estimated 
infrared luminosity. Based on analysis of a sample of lo- 
cal galaxies and assuming several galaxy SEDs, Dale et al. 
(2005) show the evolution of the ratio /tir / v /y(24^m) 
with the redshift. In our redshift range, their flgure sug- 
gests that the ratio can change by a factor of 2 - 3 for their 
galaxy sample. However, our galaxies are all luminous star- 
burst galaxies. Smith et al. (2007) estimated that the ra- 
tio of the PAH-to-TIR luminosities for starburst galaxies 
L{PAH)/Ltir ~ 0.11 ± 0.04. Marcillac et al. (2006) con- 
sider the scatter of the mid-IR indicators against each other 
and the radio (which should be a measure of the far in- 
frared) and demonstrate that, for the best fltting templates, 
the scatter is only of order 30-40% rms. Using the relative 
band strengths in Smith et al., we can estimate that only 
~ 25% of the flux density in the MIPS 24^m band is due to 
aromatic features in our redshift range. Since we apply the 
rest-frame 12pm calibration from Chary & Elbaz (2001) to 
estimate Ltir from v /y(24pm), we include the aromatic 
contribution to first order. Therefore, the relative contribu- 
tion of the aromatic features to the 24pm signal is not a 
major issue. As we show later, the ratio for our two galaxies 
detected at 70pm, /i, (70pm) //i, (24pm) « 15, is consistent 
with warm dust temperatures and of the same order as typ- 
ical starburst galaxies such as Mark 33 in Dale et al. (2005) . 
We therefore expect only a modest dispersion in the calibra- 
tion from V fv{2Afj,m) to Ltir (e.g. Papovich et al. 2002). 

More crucial is the possibility that aromatic features 
move in and out of the 24pm MIPS band in the 0.9 ^ 2 ^ 1.3 
redshift range, creating a variable calibration. The main fea- 
tures of interest are at 11.3 pm, 12.0 pm and 12.6 pm. Any 
variability due to the first and second features is negligible 
since they remain within the MIPS 24pm band over almost 
all our redshift range. However, the 12.6 pm feature (plus 
the 12.8 pm [Nell] line) moves out of the filter at z ~ 1.1. To 
check whether this produces a calibration change (namely a 
decrease of the estimated luminosity above z — 1.1), we di- 
vide our sample into two ranges below and above z — 1.1. 
The ratio of the 24pm flux densities is about unity in the 
two redshift ranges. This implies an increase of SFRtir by 
a factor of ~ 1.5 with the redshift due to geometrical effects, 
which is consistent with our calculations. The average dust 
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attenuations we calculate in the two adjacent redshift ranges 
are constant. We conclude that the derived behavior of the 
galaxies has no discontinuity at z ~ 1.1 and that there is no 
strong variation due to the aromatic feature moving out of 
the 24/im band. Therefore, we will assume that we can trust 
our estimated Ltir- 

We apply Kennicutt's (1998) calibration of Ltir into 
SFRtir, which is valid for starbursts: 

SFRtir =4.5 X W-**Ltir [erg.s~^]MQ.yr~^ .Mpc^ , 

As already noticed in Sect. 2, our observed LBG sample 
is not complete below NUV = 24.8 because of the FUV 80% 
limiting magnitude at FUV = 26.8. To estimate the total 
SFRs and surface densities down to NUV = 26.2, we assume 
that the characteristic percentages evaluated where we are 
complete can be applied below NUV — 24.8. Namely, that 
means that we assume, at NUV < 24.8 and NUV ^ 24.8, 
the same percentages of LBGs in the UV-selected sample, 
the same percentages of BLBGs and RLBGs in the LBG 
sample and identical contributions to SFRs for the above 
classes. In other words, we make the assumption that the 
instrumental threshold is not related to a physical one and 
that no important trends exist as a function of the FUV 
luminosity Lfuv- This is very likely not true over a wide 
range of Lpuv but the two Lpuv ranges (median, standard 
deviation) above and below NUV — 24.8 are similar at 1.4 x 
10^°±1.1 X lO^^Lo and 0.5x 10^°±0.2x 1O^°L0, respectively. 

Although they are more numerous, BLBGs undetected 
at 24^m do not contribute much to the SFRdust budget 
but only to SFRium one; they represent 23% of the LBG 
SFRtot. RLBGs provide most of the LBG SFRtot (77%) 
deduced from both the UV and infrared emission. However 
the surface densities of BLBGs and RLBGs follow an op- 
posite trend: 83 % arc BLBGs and 17 % arc RLBGs. Most 
LBGs are therefore blue and contribute to a small percent- 
age of the LBG SFRtot while only a small proportion of 
dusty LBGs contribute the bulk of the LBG SFRtot- 

At 1150A, LBGs furnish 69% of SFRii-,o. The ratio 
SFRtir/SFRii5o = 24.6 is larger than the ratio of SFR 
densities ptir/puv = 5.2 found at A = 1500A and in the 
redshift range 0.8 ^ ^ ^ 1.2 by Takeuchi, Buat & Bur- 
garella (2005). However, the data in Table 6 are integrated 
over the total IR luminosity range Ltir > W^'^Lq and 
iiiBO > 5 X 10^ L0. In addition to the wavelength difference, 
this might explain the different ratio. Indeed, if we sum up 
our TIR SFR density ptir = O.IIMq yr'^ Mpc'^. For 
comparison, we can estimate SFRtir. for the IR-selectcd 
sample of LIRGs and ULIRGs defined by Le Floc'h et al. 
(2005) in the same CDFS field down to SSp-Jy. The result 
is very consistent: ptir = O.I2M0 yr~^ Mpc~^ . On the 
other hand, from the FUV luminosity density provided by 
Schiminovich ct al. (2005) for Lfuv > Lmin = 0.2L*,z^3, 
we can estimate SFR1500 = 0. OO45M0 .lyr"^ Mpc~^ in 
the same FUV luminosity range, which is very consistent 
with our sum SFRnso = O.OO46M0 yr-^ Mpc'^ in the 
UV luminosity Lnso > Lmin = 0.1L*,z^3 estimated from 
NUViim = 26.2. From the SFR densities provided by Le 
Floc'h et al. (2005) and Schminovich et al. (2005), we can 
compute: Ptir/ Pfuv = 26.6, which is very close to our 
ratio. A more complete analysis requires a more sophisti- 



cated analysis using bivariate LFs that will be performed in 
a forthcoming paper (Takeuchi et al. 2007). 



6 CAN WE ESTIMATE TOTAL STAR 

FORMATION RATES OF UV-SELECTED 
GALAXIES FROM UV AND/OR IR ? 

Since the early works by Steidel et al. (1996) and Madau 
et al. (1996), astronomers have tried to evaluate the cosmic 
star formation rate and its evolution with redshift from UV- 
selected galaxies. However, a correction for the dust attenu- 
ation must be applied to the UV measurements to obtain the 
total star formation density. The most widely used method is 
based on the Meurer, Heckman & Calzetti (1999) IRX — 13 
relation, whose validity is still debated (e.g. Goldador ct 
al. 2002, Buat et al. 2005). Other methods have also been 
proposed (e.g. Kong et al. 2004, Burgarella et al. 2006a). 
Burgarella et al.'s (2006a) two-colour method, estimated for 
a sample of normal and local galaxies, does not seem to 
provide good dust attcrmation estimates, which might indi- 
cate different dust attenuation laws for these local galaxies 
than for higher redshift star-forming LBGs. Burgarella et 
al. (2006b) presented a comparison of SFRtot (i.e. FUV 
-I- dust) estimated from the FUV uncorrected for dust at- 
tenuation, the FUV corrected for dust attenuation using the 
IRX — /3 method and from IR emission. In this paper, we 
update and complete this picture (Fig. 5). The main re- 
sult is that IR-based SFRs provide the best estimate for 
SFRtot of RLBGs, which are the only LBGs for which 
we have a direct estimate of the total SFRs. This was ex- 
pected since RLBGs are dominated by their IR emission. It 
will be interesting, though, to extend the diagram to lower 
SFRs to follow the evolution. Iglcsias-Paramo ct al. (2006) 
performed a similar comparison at a ~ and found that 
SFRfuv and SFR^ust contribute about the same percent- 
age to SFRtot « 15M0 yr~^ which is not inconsistent 
with an extrapolation from our present data. On the other 
hand, the IRX — B method seems to provide, on average, 
a correct estimate of SFRtot, even for IR-bright LBGs. 
However, the dispersion is much larger that for FIR-based 
SFRs. 



7 THE EVOLUTION OF THE TOTAL IR TO 
FAR UV LUMINOSITIES OF LYMAN 
BREAK GALAXIES 

Buat ct al. (2007b) defined a sample of 190 LIRGs at 
z ~ 0.7 in the CDFS and compared it to a sample of 
120 LIRGs at 2 ~ extracted from Buat et al. (2005). 
They found that the distribution of Ltir/Lfuv is differ- 
ent at the 95% level, as a result of different average FUV 
dust attenuations: < Afuv{z = 0) >= 3.8 ± 0.1 and 
< Afuv(z ^ 0.7) >= 3.4 ± 0.1. Since the luminosity func- 
tion of LBGs is statistically consistent with the luminosity 
function of a UV selection, we can compare Ltir/Lfuv for 
our objects to the UV-selected sample at z = of Buat ct 
al. (2007a) and the BM/BX FUV-selected sample at ^ 2 
of Reddy et al. (2006). Unfortunately, Buat et al.'s sample 
does not extend far into the 10^^ regime where galaxies 
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Table 6. Contribution of 1150A UV SFR {SFR1150) and dust SFR (SFRaust) of LBGs and other non-LBG objects to the total SFR. 
SFR1150 is computed from observed GALEX NUV magnitudes (i.e. rest-frame FUV) whereas SFR^j^^t is computed from observed 
Spitzer/MIPS 2Afj,m flux density, transformed into Ldust using Chary & Elbaz (2001) and into SFRaust using Kennicutt's (1998) 
relation. 





SFRisG 
(BLBGs; RLBGs) 


SFRnoi-LBG 

(UV; IR) 


SFRall 


SFR1150 (M0.j/r-l.deg-2) 
SFRjh (M Q .yr-Kdeg-2) 
SFRtot (M .yr-^.deg-^) 
Surface Density {deg~^) 


11086 (78 %; 22 %) 
27029 (0 %; 100 %) 
38115 (23 %; 77 %) 
6003 (83 %; 17 %) 


4937; 
0; 365445 

372382 
1733;4312 


16023 
394474 
410497 
6851; 5197 




1 10 100 1000 

Figure 5. For LBGs detected at 24/xm, with Spitzer/MIPS, SFRs estimated from the dust emission are the best tracer of total SFRs. 
However, those objects are IR bright and this good correlation is expected. SFRs estimated from the UV without correcting for the dust 
attenuation underestimate the total SFRs. Finally, SFRs evaluated by applying the IRX — /3 method give a rough order of magnitude 
for the total SFRs but the dispersion is much larger than for IR-based SFRs. 



lie at z ~ 1 and z ~ 2, making a direct comparison impos- 
sible. However, Buat et al.'s FUV-selected and FIR-selected 
samples do not show any difference up to Loq^Ltot) ~ H- 
Therefore, we adopt two values at z = for the sake of 
a redshift comparison of the {Ldust/ Lfuv) ratio: first, we 
assume that we can extrapolate the FUV trend to higher 
Loq^Ltot) through a linear regression and second we use 
the FIR value. At LoqILtot) = 11-5 (i.e. a LIRG), the typ- 
ical value of Ltir/ Lfuv for the volume-corrected sample 
at 2 = is 1.92 ± 0.52 for the FIR-selected galaxies and 
1.53 ±0.36 for the FUV-selected ones. The comparative val- 
ues are 1.22 ± 0.23 at 2 = 1 and 0.84 ± 0.22 aX z = 2 (Fig. 
6). The corresponding FUV dust attenuations, using the 
Ltir/ Lfuv to Apuv calibration provided by Burgarella 
et al. (2006a), are respectively Afuv ~ 3.9, 3.0, 2.4 and 1.7. 

We clearly observe a regular decrease of Ltir / Lfuv 
with redshift (Fig. 7). What is the evolution beyond 2 = 2? 



We need more data at all luminosities and, of course, a better 
coverage of the redshift scale. The evolution of this ratio will 
be a very powerful tool to bring us some constraint on the 
dust formation timescale and process in galaxies with, later 
on, some idea of where in the UV or FIR we will have to 
look for high redshift star formation. Ultimately, dust-free 
star formation might be observed at an unknown redshift 
depending on the time needed to form dust in star forming 
regions. We must stress, however, that this decrease is based 
on a UV-selected sample (which is generally associated with 
lower stellar masses than IR-selected galaxies). Buat et al. 
(2007) find a smaller effect for an IR-selected sample at 2: « 
0.7, which might suggest a mass-dependent phenomenon. 
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Figure 6. Evolution of L^^i^t/ Lpjjv ratio for UV-selected samples as a function of the redshift from 2 ~ (Buat et al. 2005) to 2 ~ 1 
(this sample: filled dots, heavy continuous line and small crosses) and to z ^ 2 (Reddy et al. 2006: thin solid line). At 2 = 0, the linear 
regression follows the same range in Log{LxoT) covered by the FIR selection (dashed line) in Buat et al., whereas the UV selection 
(dot-dashed line) is the linear regression extrapolated beyond the observed limit at Log{Lj'OT) ~ 11.2. The dotted line indicates the 
boundaries above which we cannot observe LBGs because of the GALEX and Spitzer sensitivity limits. Even though these boundaries 
limit the Log{LTOT) range (small crosses), they should not have any effect on the J^duatl ^FUV range. 



8 THE DETECTABILITY OF LYMAN BREAK 
GALAXIES 

Will we be able to detect more UV-selected galaxies and 
LBGs at high redshift in the future? The main limitation for 
deep observations in the FIR are due to (i) thermal noise, 
(ii) cirrus and (iii) confusion limits (Helou & Beichman 1990, 
Dole et al. 2004, Kiss, Klaas & Lemke 2005). Using the in- 
formation available on the next generation of FIR/submm 
instruments, we can assess whether they would be able to 
detect LBGs. In Fig. 8, we plot the average SED of RLBGs 
and BLBGs as determined in this paper (Table 4). 

On top of our two template SEDs, the flux density lim- 
its of several telescopes are also drawn to evaluate the de- 
tectabihty of a larger sample of 2; ~ 1 — 3 LBGs/UV-selected 
galaxies in the future. The comparison is based on a Ih ex- 
posure time with a So" detection for all the instruments. The 
Atacama Large Millimeter Array (ALMA) can detect RL- 
BGs up to 2: = 2.5 in Ih and can just detect BLBGs sX z = \. 
A positive K-correction applies in the millimeter range but 
those objects are fainter even at low redshift and the K- 
correction does not help. A longer exposure time is neces- 
sary for detections at z > 1. The Herschel space telescope 
can detect the brightest RLBGs around lOO/im but will not 
be able to detect fainter ones or higher redshift ones in Ih. 
Deeper exposure times are also needed to broaden the spec- 
tral range. However, important information will be provided 
at the tip of the dust emission SED to further constrain 



the dust temperatures and TIR luminosities. The European 
SPICA Instrument (ESI) on the Japanese-led cryogenically- 
cooled telescope SPICA can take advantage of the low back- 
ground to benefit from the low confusion-limit for a 4-m FIR 
telescope. We assume, here, two kinds of detectors depend- 
ing on which one will be available for ESI in the mid-2010 
frame: Spitzer- or Herschel-like photo-detectors (GOAL) or 
state-of-the-art photodetectors (SOAP) utilizing bolometers 
with transition edge superconductor temperature sensors 
(TES) (Swinyard et al. 2006). ESI would do a fair job in 
Ih below 80/im for photoconductors and reach exceptional 
sensitivies for the TES-based devices. In the first case, it 
would provide us with individual data on z = 1 BLBGs 
while in the second case, all LBGs would be detected up to 
z — 2.0 in Ih. Finally on the blue side of the dust emission 
JWST/MIRI would be able to detect the mid-IR emission of 
low-redshift LBGs (z of 1 - 2) but higher redshift ones would 
be below the detection limits as the aromatic features shift 
out of the longest photometric band. 



9 SUMMARY 

We applied a PSF-fitting method to estimate the fiux of a 
sample of NUV=26.2 UV-selected objects in the GALEX 
Deep Imaging Survey of the Chandra Deep Field South. 
The GALEX deep images start to be confusion-limited (~ 
20 beams per source) and fluxes of objects fainter than 
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Figure 8. The expected performances of JWST, SPICA/ESI, Herschel and ALMA are overplottcd (thick Uncs with ranges of corre- 
sponding facihties quoted at the bottom of the figure) in this diagram over our template SED for BLBGs (solid lines, the higher the 
redshift the thinner the line) and for RLBG (dashed lines, the thinner the redshift the thinner the line). GOAL uses the characteristics 
of presently available photoconductors (e.g. Herschel) whereas the State Of the Art Photodetectors (SOAP) may be available in the 
mid-2010 time frame. 



NUV ~ 22 are over-estimated by the GALEX pipeline 
because objects can be blended. The selection of galaxies 
through a FUV dropout provides a very complete sample of 
star-forming galaxies in the redshift range 0.9 ^ z ^ 1.3. 
However, the efficiency is low as compared to a color-color 
method (e.g. Adelberger et al. 2004). 

We analysed the spectral energy distributions and the 
luminosity function of this sample of z ~ 1 LBGs, us- 
ing both GALEX and Spitzer/MIPS data. We divide the 
sample into two sub-classes depending on whether they are 
detected at 24/^m (red LBGs, RLBGs) or are not (blue 
LBGs, BLBGs). The NIR/MIR part of both SEDs indi- 
cates that AGNs are unlikely to play a major role in these 
objects; they are dominated by star formation. The UV 
to NIR part of the two template SEDs can be fitted by 
PEGASE models with exponentially decaying star forma- 
tion histories of ages 500 Myrs and 250 Myrs for RLBGs 
and BLBGs respectively and the same dust attenuation law 
proportional to (A/Ap-jj^))""'^- The only difference comes 
from the amount of dust attenuation: Apuv ~ 2.5, con- 
strained by the Ltir/ Lfuv for RLBGs and Afuv ~ 1.8 for 
BLBGs. The mean Ltir = lO"'^ for RLBGs (i.e. LIRGs) 
and Ltir = lO^"'^ for BLBGs. 



The BLBG SED template at z ~ 1 has a shape similar 
to the mean SED of LBGs at z ~ 3 published by Forster- 
Schreiber et al. (2004). This means that the latter objects 
undergo low dust attenuation, as expected from the selec- 
tion. However, observations by Huang et al. (2005) suggest 
that a population of dusty LBGs with large total SFRs also 
exists at high redshift. Such a population appears to be less 
numerous than the BLBGs, as is also the case at 2; ~ 1. 
On the one hand, we find that the summed UV star for- 
mation rate (uncorrected for dust attenuation) of our LBGs 
only represents ~ 1/3 of the total (i.e. UV -|- dust) star for- 
mation rate of all LBGs. That implies a large and uncertain 
extrapolation to get the total star formation density. On the 
other hand BLBGs, not detected in FIR, form the bulk (in 
number) of the galaxy population selected in UV. 

Our LBG sample at 2; ~ 1 has a UV luminosity func- 
tion consistent with that of purely UV-selected galaxies. The 
similarity of the two luminosity functions implies that mak- 
ing use of color selection is not useful to study star-forming 
galaxies whenever redshifts are available. UV measurements 
corrected using th IRX — f3 method provide correct (on 
average) but relatively uncertain total SFRs, even for our 
LBG / LIRG sample. The dispersion is larger than an es- 
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Figure 7. Prom the present data, we seem to observe a regular 
decrease of LtirI T-'FUV with redshift from 2 = (Buat et al. 
(2007a) to z = 2 (Reddy et al. (2006). Both an exponentially 
decaying law (dotted line) and a linear regression (dashed line) are 
a fair representation of the decline of the FUV dust attenuation 
in this redshift range. 



timate involving both UV and IR measurements. For RL- 
BGs, SFRtir provides the a reasonably accurate estimate 
of SFRtot. 

Finally, we observe an apparent regular decrease of the 
ratio Ltir/ Lfuv for UV-selected galaxies from 2 = to 
z ^ 2, using Buat et al. (2005) and Reddy et al. (2006) 
values in addition to ours. If confirmed, this trend might 
have strong cosmological implications in terms of the star 
formation history of the universe and the timescales for dust 
formation in primordial galaxies. 
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